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When ferromagnetic particles are suspended at an interface under magnetic fields, dipole-dipole 
interactions compete with capillary attraction. This combination of forces has recently given promis¬ 
ing results towards controllable self-assemblies, as well as low Reynolds swimming systems. The 
elementary unit of these assemblies is a pair of particles. Although equilibrium properties of this 
interaction are well described, dynamics remain unclear. In this letter, the properties of magneto¬ 
capillary bonds are determined by probing them with magnetic perturbations. Two deformation 
modes are evidenced and discussed. These modes exhibit resonances whose frequencies can be de¬ 
tuned to generate non-reciprocal motion. A model is proposed which can become the basis for 
elaborate collective behaviours. 
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I. INTRODUCTION 

The deformation of a liquid surface by identical float¬ 
ing particles [1( makes them attract each other mm- a 
practical way to avoid clustering consists in using parti¬ 
cles possessing a large magnetic moment m HE], if m is 
perpendicular to the interface, the dipole-dipole interac¬ 
tion is repulsive [6j and opposes the capillary attraction. 
The combination of a repulsive magnetic interaction with 
the attractive capillary interaction creates a pair poten¬ 
tial possessing a minimum. Two particles will then settle 
at an equilibrium distance creating a ‘magnetocapillary 
bond’ which is the building block for larger stable struc¬ 
tures HIZ]. Indeed the energy scale associated with this 
magnetocapillary interaction can be much larger than the 
thermal energy for a wide range of submillimetre-sized 
particles. Particle size is bounded from below by the 
capillary attraction, which decreases as particles become 
smaller. At room temperature, the size at which ther¬ 
mal energy becomes strong enough to break free from 
the capillary interaction is 3.4 fim in diameter, as can be 
estimated from the potential given in equation [l] How¬ 
ever, one could generate a stronger capillary attraction 
by using confined geometries mm or other body forces 
than gravity pQ. To allow remote tuning of the bond 
properties, it has been proposed to use particles with 
large magnetic susceptibility x i n an external magnetic 
induction B uni. In this case fh = xVBlfi o (with V the 
volume of the bead), meaning that varying the magnetic 
induction allows to control the strength and direction of 
the magnetocapillary interaction [5:. Adding an oscillat¬ 
ing horizontal field deforms these aggregates [10 which 
are observed to swim at low Reynolds number m by 
spontaneously changing their shape in a non-reciprocal 
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FIG. 1. Pair potential. - a) Notations. Two beads are 
placed at a water surface and submitted to a constant induc¬ 
tion B — B x e x + B z e z as well as a small oscillating induction 
p. Distance d and angle 0 characterise the system, b) Mag¬ 
netocapillary potential. Interaction potential U — U c -\-Um 
in polar coordinates ( d , 0 ), as given by eq. for typical pa¬ 
rameters of the experiment. Darker zones indicate potential 
wells. On the left half, vertical induction B z generates an 
isotropic repulsion when Me s= 0, which leads to a finite equi¬ 
librium distance d eq ~ 2.5 mm. On the right half, cylindrical 
symmetry is broken when Me / 0. Me is the orthoradial mag¬ 
netocapillary number introduced in eq. ©• It is proportional 
to Bf.. 


way and, what is more, swimming direction can be re¬ 
motely controlled [12;. 

The aim of the present letter is to study experimen¬ 
tally the fundamentals of the magnetocapillary bonds 
and their vibration modes. First, we present their poten¬ 
tial energy and their static equilibrium distances. Then, 
their two normal vibrations (stretching and swinging) are 
evidenced experimentally as well as modelled in a pertur¬ 
bative way. Finally, we show that the detuning of their 
resonance frequencies can be exploited for creating non¬ 
reciprocal cycles at low Reynolds number. Although par- 
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tides of a few hundred microns in diameter were used, 
if the system is downscaled, resonance of the swinging 
mode can still be achieved. This means that similar dy¬ 
namics could be reproduced with smaller, colloidal par¬ 
ticles [13] . 


around 0 = 0 and the equilibrium distance d eq . It can be 
shown that d eq = d eq /l c depends only on Md following 
the equation 

d* g K 1 (d eq ) = 3M d . (4) 


II. MAGNETOCAPILLARY PAIR POTENTIAL 

We consider two floating beads of diameter D possess¬ 
ing identical magnetic moments fh = mn x e x + m z e* zi e x 
and e z being the horizontal and vertical unitary vectors, 
respectively, as defined in fig. [l]a) . Their center-to-center 
vector is denoted d = de r , and the orientation of the pair 
relatively to e x is 0. Their interaction can be modelled by 
a pair potential U composed of the sum of an attractive 
term due to capillarity U c and a magnetic dipole-dipole 
U m term that can be either attractive or repulsive de¬ 
pending on HDHH]. 

A reasoning similar to [ 2 ] will be used to express the 
capillary interaction energy U c . Solving Laplace equa¬ 
tion for a single sphere gives interfacial shape ((d) oc 
Kq ( d/l c ), where l c = y/j/pg is the capillary length and 
Kq is the modified Bessel function of the second kind 
and order 0. If d D, the superposition approxima¬ 
tion can be used to express U c , meaning that the total 
interfacial deformation is the sum of the independant de¬ 
formations caused by each particle [14]. As will be shown 
in the discussion of figure [4] this approximation is in good 
agreement with the experiment for d > 2D. We find 


If Mq = 0, the potential is axisymmetric and every ori¬ 
entation 6 is equivalent for the pair. Figure [l] b) shows 
the interaction potential U in polar coordinates, with one 
bead at the axis origin, for typical values of the experi¬ 
ment. Because the potential is mirror symmetric around 
the 90° axis, the two halves of the graph are chosen to 
represent two cases. Case Mq = 0 is on the left, with 
a minimum of energy at a given distance d eq . On the 
right, the addition of a horizontal induction B x of order 
10 -3 T ( Mq 7 ^ 0) breaks axisymmetry and creates an 
equilibrium orientation at 0 = 0. Typical values for the 
experiment are B z « 5 x 10 -3 T and B x « 2 x 10 -3 T, 
which gives a ratio Mq/Mq « 0.8. 

Using a perturbative analysis of [/, the stiffness of the 
potential well can be obtained for its two principal direc¬ 
tion (radial kd and orthoradial ke): 


kd 


r 

rfd 


6Ki(d) - 2 dK 0 (d) 


( 5 ) 


kg 



l%d 5 


(6) 


U c = -TK 0 (d/l c ) (1) 

with r = 2ir / yq 2 being a typical capillary energy scale. 
The capillary charge of a particle q is a length that char¬ 
acterises the interface deformation needed to compensate 
the buoyancy of the particle. It depends explicitly on 
the bead volume, the density of the bead, the density of 
the liquid, the capillary length and the wetting contact 
angle [2 , 3 . If none of the fluids has strong magnetic 
properties, the magnetic dipole-dipole potential is 

tt _ no [ra 2 + ra 2 (l -3 cos 2 6>)] , oA 

Um ~ 4t r# W 

with /iq the void permeability. If the distances are nondi- 
mensionalised by the capillary length (d = d/l c ), the 
shape of pair potential U = U c + Um is determined by 
two dimensionless magneto capillary numbers Md and Mq 
such that 


u = r 


TT _L Md _L • 2/j 

—Ko (d) + sin 6 

d 6 d 6 


( 3 ) 


where Md = /io — 2 m^)/ 47 rT / 3 is the radial magneto¬ 
capillary number and Mq = 3/io^/47rT/ 3 the orthora¬ 
dial magneto capillary number. 

While the magnetic terms of U are dominating at both 
short and long distances, it has a stable minimum pro¬ 
vided that Md < 1.11. Thus, it forms a potential well 


In the case of two identical beads of mass W\ = W 2 = 
W, the reduced mass is W/2. Two frequencies can thus 
be associated with the magnetocapillary well: a radial 
frequency 


Ud = Znfd = ^/ZhJW 


2 Ucideg) 9M d r 

lc\ W + Wd%’ 


( 7 ) 


assuming d^> D, that corresponds to a stretching mode 
and an orthoradial frequency 


LOg = 2nfg = \f2dzgJW 


2 Mg I 2r 

TT ^Wd% q 


(8) 


that corresponds to a swinging mode. At low Reynolds, 
the attenuation is due to the viscosity of the ambient 
fluid, so that attenuation time is 


W W 
~ r]D 


( 9 ) 


where C is the damping coefficient. We have Cstokes = 
37 tt]D in the case of a fully immersed sphere mansi. if 
fdJ > 1 or for > 1, we can force harmonic resonances, 
even in the Stokes regime. This is of interest as frequency 
can thus be changed to adjust phase. 
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FIG. 2. Characterisation of the system. - a) Beads 
magnetisation. Magnetic field H induces a magnetisation 
m in a 500 /im bead with negligible hysteretic behaviour. H is 
successively increased from 0 (1), decreased (2) and increased 
back to 0 (3). Values of H are typical of the experiment. In¬ 
set shows a larger cycle up to saturation magnetisation, with 
same units, b) Equilibrium distance. The dimensionless 
equilibrium distance d eq as a function of the radial magne¬ 
tocapillary number Md . The colourmap shows distance d eq 
over diameter D. Solid line corresponds to theoretical predic¬ 
tion without adjustable parameter, obtained by numerically 
solving eq. Symbols represent different bead diameters. 


III. EXPERIMENTAL SETUP 

In order to study the vibrations of magnetocapillary 
bonds, we will use the following experimental setup: a 
large Petri dish is filled with water and placed at the cen¬ 
ter of a three-axis, Earth’s field compensating Helmholtz 
coil system. Two identical submillimetric chrome steel 
beads (diameter D = 397, 500 or 793 /im, density p = 
7830 kg/m 3 , relative magnetic susceptibility x > 300) 
are disposed at the water-air interface. The capillary 
charges of the beads, calculated for a contact angle of 
90°, are 6, 12 and 45 /im, for the respective diameters 
397, 500 and 793 pm. In the presence of a constant mag¬ 
netic induction B = B x e x + B z e z , the beads settle at an 
equilibrium distance d eq . For vibration experiments, a 
small oscillating induction (3 = /3(t) (cos sin aej) is 

added. Angle a can be varied to produce different oscil¬ 
lating behaviours. The kinematics of the system is fully 
described by the beads center-to-center distance d and 
its rotation angle 0 as defined in fig. [l] a). 

For obtaining the desired magnetic properties, we 
chose chrome steel spheres (alloy AISI 52100). The bulk 
material is ferromagnetic and has a large susceptibility 
X, which depends on applied magnetic field H = B/pq, 
such that x(iJ) > 300 [17]. However, a finite ferro¬ 
magnetic body in a magnetic field produces a demag¬ 
netisation effect. An effective susceptibility Xeff takes 
this effect into account, which for a sphere has the form 
Xeff = x/(l + x/3) m For X > 3, the sphere will 
thus behave like an isotropic superparamagnetic particle 
of susceptibility Xeft = 3. In fig. [2] a), we plot the experi¬ 
mental magnetisation curve m(H). Indeed, no hysteresis 
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FIG. 3. Resonance spectrum. - a) Amplitude. Fre¬ 
quency response of distance d and angle 0 for typical values of 
the experimental parameters (D — 500 pm , Me/Md ~ 0.8). 
Radial oscillation of amplitude Ad happens for a = 0°, while 
orthoradial oscillation of amplitude A0 happens for a — 90°. 
Error bars are smaller than the symbols, with a typical error 
of 0.2 pm on the bead center. Data have been fitted with a 
resonance curve, represented by solid lines, b) Phase. The¬ 
oretical prediction for the phase in both cases. Colourmap 
indicates phase difference A 


is observed, and the slope gives XeS — 3 as expected. A 
small magnetisation remains at zero field, of the order of 
100 Am -1 , which is negligible for the fields used herein 
(typically 5000 Am -1 ). 


IV. RESULTS 

Equilibrium distance d eq has been measured for 3 di¬ 
ameters D and various values of radial magnetocapillary 
number M^, changed by varying both B x and B z . As 
seen in fig. |2]b) , the equilibrium distance calculated from 
the interaction potential, as defined in eq. 0, is in good 
agreement with experimental values for large distance d eq 
compared to diameter D. Indeed, for d eq /D < 2 (light- 
coloured points in the graph), the superposition approx¬ 
imation no longer holds in the determination of capillary 
potential U c B3|. In this case, the attraction is stronger. 

In order to characterise the properties of the magne¬ 
tocapillary well, a horizontal oscillating induction /3(t) 
with magnitude /3(t) = /3sin(27r/t) and angle a is added. 
A small oscillation amplitude (3 <C \B X \ is chosen, so 
that the system remains close to its equilibrium configu¬ 
ration. Pertubative analysis gives two oscillation modes 
in radial and orthoradial directions, that we will excite 
independently. Two values of angle a corresponding to 
both directions will thus be studied: a = 0° (radial) and 
a = 90° (orthoradial). 

In fig. [3] a), oscillation amplitudes of distance d and 
angle 0 are represented as a function of frequency. Each 
point was obtained by identifying the geometric center of 
each bead, averaged over 10 oscillation periods, in order 
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FIG. 4. Resonance frequencies. - Comparison between 
measured resonance frequency f exp and theoretical prediction 
fth , for several values of induction fields B x , B z and /3 as 
well as several bead diameters. Two orientations a of fd are 
explored which correspond to radial and angular oscillations, 
respectively a = 0° on the left and a = 90° on the right. The¬ 
ory and experiment agree well, with the exception of distances 
d/D < 2 in case a = 0°, corresponding to light-coloured 
points in the left graph. 
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FIG. 5. Deformation cycles. - Trajectories in plane 
(0, d/D) for 10 periods of the oscillating induction with a m 
45°. Below both resonance frequencies, in-phase oscillations 
produce no loop. When the frequency is increased, loops ap¬ 
pear as oscillations get out of phase. Above both resonances, 
cycles disappear as phase difference decreases to zero. 


to reach sub-pixel precision. Two independent oscillation 
modes can indeed be excited: a radial mode of amplitude 
Ad for a = 0°, for which d oscillates around d eq and 
0 = 0 (stretching mode), and an orthoradial mode of 
amplitude A0 for a = 90°, called swinging mode, for 
which 0 oscillates around 0 and d = d eq . Both modes 
exhibit a resonance, with distinct resonance frequencies 
fd > fe in this example. 

Steady-state solution of a forced damped harmonic 
oscillator is usually defined as the real part of a com¬ 
plex number z(t) = R(uj)e luJt with function R(uj) = 
If Fo/W is the driving force per unit mass, 
we have the following expressions for amplitude A(ui) and 
phase ip(uj), respectively 


A(lu) — 


Fp/W 

^/ ((jJq — ^ 2 ) 2 + 4c0> 2 /r 2 


( 10 ) 


and 


V , M=«rctan( T( ^ ij2) ). (11) 

with angular frequency u) = 27r/, resonance frequency 
ujq = 27t/o and damping time r. Data in fig. [ 3 ] a) 
are very well fitted by function A(uj), with fit parame¬ 
ters Fq/VF, cjo and r. The resonance curves correspond 
to underdamped harmonic oscillators as the quality fac¬ 
tor Q = ujot/2 ranged between 0.6 and 3.8. The reso¬ 
nance frequencies in all our experiments ranged between 
0.5 and 5.2 Hz. Finally, the time scale r given by the 
fit corresponds to an average viscous damping coeffi¬ 
cient C that grows linearly with bead diameter D and 
is close to a Stokes damping Cstokes = 37 rr]D, such that 
C/Cstokes = 0.86 zb 0.04. This seems consistent with the 
viscous damping of a partially immersed sphere mm- 


The dispersion of the values can be due to the interac¬ 
tion between the spheres, as well as slight variations of 
the wetting conditions. Note that this average value of C 
excludes radial oscillations for which d/D < 2, which are 
discussed below. Figure [ 3 ] b) represents the correspond¬ 
ing phases ip. Because radial and orthoradial resonance 
frequencies are distinct, a phase difference Aip can arise 
between modes. Close to resonance frequencies, we have 
Aip> ~ 7 t / 2 . 

Radial and orthoradial resonance frequencies, respec¬ 
tively fd and /#, are determined for different sets of pa¬ 
rameters B x , B z , /3 and D. Figure [4] compares experi¬ 
mental values fd , exp an d fe, exp with theoretical predic¬ 
tions fd , th an d fe, th as defined in eqs. (J7| and (J8|. Case 
a = 0° shows good agreement for large values of d/D , 
but theory gives an overestimation of frequency fd for 
distances d/D < 2. This evidences the limits of vali¬ 
dation of the superposition approximation in capillary 
potential U cap . In case a = 90°, by contrast, theory and 
experiment agree well. Indeed, linear regression gives 
fd, exp ~ fd, th down to about 1%. Because distance d 
is constant in the orthoradial case, capillarity does not 
act as a spring force, which explains why this limit of 
validation is not observed. 


V. DISCUSSION 

Controlling the phase difference between degrees of 
freedom is the key ingredient for generating significant 
non-reciprocal deformations, which are necessary to low 
Reynolds propulsion [2Q]. Indeed, to propel itself, a body 
immersed in a fluid must produce a net flow over a time 
period. This means that low Reynolds swimmers must 
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find a way around the time-reversal properties of Stokes 
flows, either by using the medium to break symmetry m 
or by undergoing non-reciprocal, i.e. non-time-reversible, 
deformations [1TJ 20]. Such deformations can be ob¬ 
tained in a magnetocapillary assembly by applying time- 
dependent magnetic induction fields m- Note that the 
field itself can be time-reversible, as the system sponta¬ 
neously breaks time-reversal symmetry through magne¬ 
tocapillary interactions [T2] . 

In this work, we evidenced that for magnetocapillary 
bonds, even if the fluid is in a Stokes regime, the Reynolds 
number is not the only parameter driving the dynamics 
of this system, the system inertia has to be taken into 
account. Indeed, quality factor Q = ujqt/2 can be over 
unity, meaning that the system is underdamped, granted 
that the stiffness of the potential and the particle mass 
are high enough. This is of interest for creating phase dif¬ 
ferences between stretching and swinging of the bonds. 
Each mode possesses a quality factor Q which gives the 
shape of the resonance curve, and the ratio of the fre¬ 
quency to the resonance frequency £ = uj/ujq allows to 
distinguish between the well known regimes of a forced 
harmonic oscillator. Firstly, assuming Q > 1 (for Q < 1, 
the same discussion can be made using l/Q instead of 
Q), if £ < l/Q the oscillator follows the forcing in a 
quasistatic way, in phase. Then, if Q > £ > l/Q, the 
dynamics is limited by the viscous dissipation and the 
oscillation of the system is in quadrature with the forc¬ 
ing. Finally, if £ > Q the oscillator is limited by inertia 
and oscillates in phase opposition with the forcing. 

The stretching and swinging modes of the magnetocap¬ 
illary bond are generally detuned (fd j- fo)- They can 
be both excited simultaneously using intermediate angles 
a (0 and 90° produce pure radial and orthoradial oscil¬ 
lations respectively). If the resonances are sufficiently 
separated i.e, fd/fo < l/Q or fd/fe > Q, the excitation 
frequency can be chosen so that one mode oscillates in 
quadrature and the other oscillates in phase or phase op¬ 
position with the excitation as was shown in fig. [3] There¬ 
fore, the system can follow simple non-reciprocal Lis- 
sajous curves (ellipses) in the plane (#, d/D ) whose shape 
depends on control parameters / and a. Figure [^repre¬ 
sents six cycles for a = 45°. For each cycle, 10 periods of 
the oscillating induction are represented. For frequencies 
below both resonances, i.e. f < fd and f < fo, in-phase 
oscillations produce reciprocal deformations. When fre¬ 
quency is increased, cycles appear as oscillations get out 
of phase. Around resonance frequencies, phase difference 
is close to 7r/2. At frequencies higher than both reso¬ 
nances, i.e. f > fd and f > fe, cycles disappear as 
phase difference decreases to zero. 

Further downscaling of the system could prove useful 
for future applications, like low Reynolds swimming. To 
determine if a phase difference is achievable at smaller 
scale, we can look at the expressions of quality factors 
in radial and orthoradial directions. For a constant ratio 


D/d , we have 


Qe = ~ DB - ( 12 ) 

This means that orthoradial q factor can be kept con¬ 
stant, granted that induction field B is increased accord¬ 
ingly. As a matter of fact, at the center of a Helmholtz 
system, we have B ~ 1/R where R is the radius of the 
coils. If coil size is decreased along with particle diame¬ 
ter D , with everything else kept constant, the orthoradial 
oscillator will behave similarly. In radial direction, how¬ 
ever, we have 


Qd 


oo d T qD 2 D 5 
__ - — 


(13) 


since capillary charge q ~ D 3 . With a constant ratio 
D/d, we have Qd ~ D A meaning that radial oscillation is 
overdamped at smaller scales. As a result, for small scale 
systems, the only way to achieve non-reciprocal cycles is 
to use the orthoradial resonance frequency as excitation 
frequency. 

Of course, exciting radial and orthoradial modes us¬ 
ing independent, out of phase magnetic induction fields 
would also generate non-reciprocal motion. However, 
magnetocapillary bonds in assemblies of more than 2 par¬ 
ticles are usually not aligned, meaning that both modes 
are always excited. This may explain why a single os¬ 
cillating sinusoidal field can generate efficient locomo¬ 
tion [nj HU • Previous theoretical work on magnetocap¬ 
illary systems m studied the Stokesian motion of parti¬ 
cles in a magnetocapillary well, finding insufficient sym¬ 
metry breaking to account for efficient low Reynolds loco¬ 
motion in a three beads system. However, the frequency 
was not a relevant parameter in their work. Lastly, it 
is obvious from eq. © that for larger deformations, the 
potential becomes anharmonic. Stretching and swinging 
degrees of freedom become coupled, opening a rich vari¬ 
ety of nonlinear behaviour. 


VI. CONCLUSION 

We demonstrated that magnetocapillary bonds can be 
created along liquid interfaces. Magnetocapillary num¬ 
bers Md and Mq determine the interdistance and orien¬ 
tation of the bond, respectively. Then, we proved the¬ 
oretically and experimentally that such bonds possess 
two vibration modes, a radial stretching mode and an 
orthoradial swinging mode, that can be excited indepen¬ 
dently. These modes exhibit a resonance when excited, 
even though the system is in Stokes regime. This means 
that bead inertia must be taken into account and that 
frequency is an essential control parameter of the bond. 
Indeed, depending on which frequency is applied to the 
system, it can experience out of phase radial and ortho¬ 
radial deformations. This allows to generate a breaking 
of time reversibility, which is required for low Reynolds 
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locomotion. If this system were to be downscaled, sim¬ 
ilar dynamics would occur at the orthoradial resonance 
frequency. 

While this paper focused on near equilibrium dynam¬ 
ics, nonlinear behaviours and coupling certainly matter 
for low Reynolds propulsion. Studying the response of 
magnetocapillary bonds to finite perturbations would be 
the next logical step. Further work should also deter¬ 
mine how to combine several pairs of particles in larger 
systems. Indeed, magnetocapillary bonds constitute ver¬ 
satile building blocks that can be combined to create 
elaborated microstructures. 
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